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We calculate the viscous pressure tensor of the quark-gluon plasma in strong magnetic field. 
It is azimuthally anisotropic and is characterized by five shear viscosity coefficients, four of 
which vanish when the field strength eB is much larger than the plasma temperature squared. 
We argue, that the azimuthally anisotropic viscous pressure tensor generates the transverse 
flow with asymmetry as large as 1/3, even not taking into account the collision geometry. 
We conclude, that the magnitude of the shear viscosity extracted from the experimental data 
ignoring the magnetic field must be underestimated. 

I. INTRODUCTION 

Strong magnetic field produced in relativistic heavy-ion collisions [H [2] has a strong impact 
on phenomenology of the quark-gluon plasma (QGP). It induces energy loss by fast quarks and 
charged leptons via the synchrotron radiation [3] and polarization of the fermion spectra [3]. It 
contributes to the enhancement of the dilepton production [1] and azimuthal anisotropy of the 
quark-gluon plasma (QGP) [5]. It causes dissociation of the bound states, particularly charmonia, 
via ionization (6j [7j. Additionally, the magnetic field drives the Chiral Magnetic Effect (CME) 
[U 181411). which is the generation of an electric field parallel to the magnetic one via the axial 
anomaly in the hot nuclear matter. 

It has been argued recently in [5] that the magnetic field of strength eB ~ m^. [TJ [2] is able to 
induce the azimuthal anisotropy of the order of 30% on produced particles. This conclusion was 
reached by utilizing the solution of the magneto-hydrodynamic equations in weak magnetic field. 
In this paper we discuss the magneto-hydrodynamics of the QGP in the limit of strong magnetic 
field. Our goal is to calculate the effect of the magnetic field on viscosity of the plasma. It is well- 
known that the viscous pressure tensor of magnetoactive plasma is characterized by seven viscosity 
coefficients, among which five are shear viscosities and two are bulk ones. Generally, calculation 
of the viscosities requires knowledge of the strong interaction dynamics of the QGP components. 
However, in strong enough magnetic field these interactions can be considered as a perturbation and 
viscosities can be analytically calculated using the kinetic equation. Application of this approach 
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to the non-relativistic electro-magnetic plasma is discussed in |12j . A general relativistic approach 
was developed in [IB] . We apply it in Sec. [il] to derive the viscosity coefficients of QGP, which are 



given by (22) and (33). As in the non-relativistic case, we found that four viscosities vanish as the 
magnetic field strength increases. 

A characteristic feature of the viscous pressure tensor in magnetic field is its azimuthal 
anisotropy. This anisotropy is the result of suppression of the momentum transfer in QGP in 
the direction perpendicular to the magnetic field. Its macroscopic manifestation is decrease of 
the viscous pressure tensor components in the plane perpendicular to the magnetic field, which 
coincides with the "reaction plane" in the heavy-ion phenomenology. Since Lorentz force vanishes 
in the direction parallel to the field, viscosity along that direction is not affected at all. In fact, 
the viscous pressure tensor component in the reaction plane is twice as small as the one in the 
field direction. As the result, transverse flow of QGP develops azimuthal anisotropy in presence of 
the magnetic field. Clearly, this anisotropy is completely different from the one generated by the 
anisotropic pressure gradients and exists even if the later are absent. 



In Sec. Ill we discuss QGP transverse flow in strong magnetic field using the Navie-Stokes equa- 



tions. At later times after the heavy- ion collision, flow velocity is proportional to r/ 1 > 2 , see (40a) 



and ( 40b ) . If the system is such that in absence of the magnetic field it were azimuthally symmet- 



ric, then the magnetic field induces azimuthal asymmetry of 1/3, see (44). This is surprisingly close 
to the weak field limit recently reported in [5]. The effect of the magnetic field on flow is strong 
and must be taken into account in phenomenological applications. Neglect of the contribution by 
the magnetic field leads to underestimation of the phenomenological value of viscosity extracted 
from the data |14H16j . In other words, more viscous QGP in magnetic field produces the same 
azimuthal anisotropy as a less viscous QGP in vacuum. 



II. VISCOUS PRESSURE IN STRONG MAGNETIC FIELD 
A. Kinetic equation 

Kinetic equation for the distribution function / of a quark flavor of charge ze is 

p^f = zeB^^u u + C[f,...] (1) 

where C is the collision integral and B^ v is the electro-magnetic tensor, which contains only mag- 
netic field components in the laboratory frame. Ellipsis in the argument of C indicates the distri- 
bution functions of other quark flavors and gluons (we will omit them below). The equilibrium 
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distribution: 

f n = E e -Pp-U(x) (2) 

J0 4irm 3 TK 2 (f3m) { ' 

where U{x) is the macroscopic velocity of fluid, p^ 1 = mu^ is particle momentum, (3 = 1/T and p 
is the mass density. Since oc u„, the first term on the r.h.s. of (111 vanishes in equilibrium as 



well as the collision integral. Therefore, we can write the kinetic equation as an equation for 5f 

P^d,h = zeB^ d ^u v + CW] (3) 
where 5f is a deviation from equilibrium. Differentiating Q we find 

dJo = -fo^p X d IM U x (x) (4) 
Since U x = (7y,7yY) and p x = (e,p) = (7 v m,7 v mY) it follows 



Thus, in the comoving frame 



Substituting Q in ^ yields 



where we defined 



dMv=o = f°^P»W (6) 



V^ = i(S M K + 9^) (8) 



and used u^u v d^V v = u^u^V^. 

Since the time-derivative of /o is irrelevant for the calculation of the viscosity we will drop it from 
the kinetic equation. All indexes thus become the usual three-vector ones. To avoid confusion we 
will label them by Greek letters from the beginning of the alphabet. Introducing b a p = B~ 1 E a p~ f B~ / 
we cast ([7]) in the form 

\v a u^h = -zeBb a ^^- £ - C[6f] . (9) 

The viscous pressure generated by a deviation from equilibrium is given by the tensor 

f d 3 p 

- Ha/3 = / p a pp 5f — (10) 
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Effectively it can be parameterized in terms of the viscosity coefficients as follows (we neglect bulk 
viscosities) 

1 



(11) 



n=0 



(n) 

where the linearly independent tensors V^J are given by 



a/3 



v $? = ( 3b ch - M ( WsVjs - • v 



(12a) 



V$! = 2V a p + 5 a pVy 5 b 7 b 6 - 2V a7 b 7 bp - 2V M b 1 b a + {b a bp - <5 Q/3 )V • V + b a b V lS b 7 b s (12b) 



*/3 

r(2) 



V a/3 = 2 (^a 7 &/3 7 + Vp 7 b ai - V^b^b^bs) 



(3) 

V aj3 = Vajbfij + Va 7 6a 7 — V-ysbasba'ybpbs — Vyfibp-fbabs 



V af = ZiVySbaSbaybpbs + Vysb/3jb a b 5 ) . 

For calculation of shear viscosities rj n , n = 1, . . . , 4 we can set V • Y = and V a pb a bp = 0. 

in) 

Let us expand Sf to the second order in velocities in terms of the tensors V^p as follows 
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(12c) 
(12d) 
(12e) 



(n) a 3 



(13) 



n=0 



Then, substituting (13) into (11) and requiring consistency of (10) and (11) yields 

2 



Vn 
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ev g n d p 



(14) 



This gives the viscosities in the magnetic field in terms of deviation of the distribution function 



from equilibrium. Transition to the non-relativistic limit in (14) is achieved by the replacement 
e m [T2"] . 



B. Viscosity of collisionless plasma 



In strong magnetic field we can determine g n by the method of consecutive approximations. 
Writing 5f = 5f^ + 5f^ and substituting into aty we find 



(15) 



Here we assumed that the deviation from equilibrium due to the strong magnetic field is much 
larger than due to particle collisions. The explicit form of C is determined by the strong interaction 



5 



dynamics but drops off the equation in the ieading oder. The first correction to the equilibrium 
distribution obeys the equation 



1 V f d; dS f (1)l 
^PaVpV a pf = ~ZeBb a pVp — - 



(16) 



Using (13) we get 



dSfW 

b a /3V(S-^- = 2b a pVp ^ 9n 
" n=0 



(17) 



Substituting (17) into (16) and using (12) yields: 



TzeB 



PaVpVapfo = -2bp u v a v v [gi(2V a p - 2Vp y b 1 b a ) + 2g 2 Vp 1 b 1 b a 

+gz{V on bp 1 + Vp 7 b ai - Vysb a b s ) + 2g±V~ / sbp~ f b a b s )} (18) 



where we used the following identities b a pb a = b a pbp = b a pv a vp = 0. Clearly, (18) is satisfied only 
if gi = <?2 = 0. Concerning the other two coefficients, we use the identities 

b a pbp^ = b 7 b a — 5 ai b 2 , (19a) 

£a/3 7 £<5eC = ^aS {^pe^C, ~ ^PC^e) ~ <5«e (SpgS^ - Sp^g) + <5q< {SpS^e ~ Sp^s) (19b) 



that we substitute into ( 18 ) to derive 



2TzeB 
Since p a = ev a we obtain 



p a v p V a pf = gz[2V a pb a bp - W a pv a bp{b ■ v)] + 2g A V a pv a bp{h ■ 



53 



94 e 2 /o 



(20) 



(21) 



2 ATzeB 

Using (|2j), ^ in Q in the comoving frame (of course %'s do not depend on the frame choice) 
and integrating using 3.547.9 of [21] we get 

K 3 ((3m) P T 



V3 K 2 {pm)2zeB 
The non-relativistic limit corresponds to m > T in which case we get 

„NR_ PT 

V3 ~2z~eB- 

In the opposite ultra-relativistic case m <C T (high-temperature plasma) 



zeB 



(22) 



(23) 



(24) 



where n = p/m is the number density. 
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C. Contribution of collisions 

In the relaxation-time approximation we can write the collision integral as 

C[6f) = -vSf (25) 

where v is an effective collision rate. Strong field limit means that 

ujb^v (26) 

where ojb = zeB je is the synchrotron frequency. Whether v itself is function of the field depends 
on the relation between the Larmor radius rs = vt/u>b, where vt is the particle velocity in the 
plane orthogonal to B and the Debye radius rfj. If 

r B > r D (27) 



then the effect of the field on the collision rate v can be neglected [12 . Assuming that (27) is 
satisfied the collision rate reads 

v = nva t (28) 

where at is the transport cross section, which is a function of the saturation momentum Q s |19U20| . 
We estimate at ~ a^/Q 2 , with Q s ~ 1 GeV and n = P/T with pressure a 2 s P ~ 1 GeV/fm 3 we get 
v ~ 40 MeV. Inequality ( 26 ) is well satisfied since eB ~ ml [HE] and m is in the range between the 



current and the constituent quark masses. On the other hand, applicability of the condition (27) is 



marginal and is very sensitive to the interaction details. In this section we assume that (27) holds 
in order to obtain the analytic solution. Additionally, the general condition for the applicability of 
the hydrodynamic approach I = \ jv <C L, where £ is the mean free path and L is the plasma size 
is assumed to hold. Altogether we have <C rs <C t <C L. 



Equation for the second correction to the equilibrium distribution Sf^ follows from ( 15 ) after 



substitution (25) 



Now, plugging 

= [93V^+g^f}v a v p , (30a) 
5f {2) = [9iVS + 92V$]v a vp (30b) 
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into ([29]) yields 

2zeB 

— ^— {giftVpabaryVpv-, - 2Vp a b aj v l3 Vy(Y ■ b)] + 2g 2 V l3a b ai vpv 1 {Y ■ b)} 

= -z^53 {-2^6^^^ - 6T^ a 6 Q:7 w i a'w 7 (y • b)} (31) 

where we used g$ = 2g%. It follows that 

91 = T = aJ7 (32) 

With the help of @,([2]),([T4]) we obtain 

= r/ 2 _ 8 p 2 a t T z l 2 K 7/2 (f3m) 

III. TRANSVERSE FLOW 

To illustrate the effect of the magnetic field on the viscous flow of the electrically charged 
component of the quark-gluon plasma we will assume that the flow is non-relativistic and use the 
Navie-Stokes equations that read 

/dVg | v dv a \ = bp | dii a p 

\ dt 13 dxp J dx a dx/3 
where Tl a B is the viscous pressure tensor, p = mn is mass-density and P is pressure. We will 
additionally assume that the flow is non-turbulent and that the plasma is non-compressible. The 
former assumption amounts to dropping the non-linear in velocity terms, while the later implies 
vanishing divergence of velocity 

V • V = (35) 

Because of the approximate boost invariance of the heavy-ion collisions, we can restrict our atten- 
tion to the two dimensional flow in the xz plane corresponding to the central rapidity region. 
The viscous pressure tensor in vanishing magnetic field is isotropic in the xz-plane and is given 

by 

o (dV a dV p \ I Vxx V xz \ 
where the superscript indicates absence of the magnetic field. In the opposite case of very strong 



magnetic field the viscous pressure tensor has a different form (11). Neglecting all r\ n with n > 1 
we can write 

( -Vzz \ (\V XX 

= Vol )= 277o I (37) 

V 2V ZZ / V v zz 



s 



where we also used ( 35 ) . Notice that II 



oo 

XX 



-U°° - 

2 zz 2 xx 



hll xx indicating that the plasma flows in the 



direction perpendicular to the magnetic field with twice as small viscosity as in the direction of 
the field. The later is not affected by the field at all, because the Lorentz force vanishes in the 



field direction. Substituting ( |37[ ) into ( |34[ ) we derive the following two equations characterizing the 



plasma velocity in strong magnetic field 

dV x 



dP d 2 V x 
dx ^° dx 2 ' 



dV z 



r dt dx 1 ,u dx 2 ' r dt 

Additionally we need to set two initial conditions 



V x \ t=0 = <Pl(x,z) , 



V, 



t=o 



dP 

dz dz 2 



(f 2 {x,z) 



The solution to the the problem (38), (39) is 



/OO f't /'OO 

dx'ipi{x\z)Gi{x - x',t) - - / dt' dx'Gi(x-x',t-t 
-oo 2 P Jo J-oo 2 



V z (x,z,t)= / dz'(p 2 {x,z')Gi(z - z',t) - 

J — OO 

Here the Green's function is given by 

G k (z,t) 

and the diffusion coefficient by 



P Jo 



dt' 



1 



^dP(x',z,t') 
dx' 



j i n I / i if \ dP(x, Z it) 

dz Gi(z - z ,t-t ) — 



ira 2 kt 



2r]o 
P 



(38) 

(39) 

(40a) 
(40b) 

(41) 

(42) 



Suppose that the pressure is isotropic, i.e. it depends on the coordinates x,z only via the radial 



coordinate r = \J x 2 + z 2 ; accordingly we pass from the integration variables x' and z' to r in (40a) 



and (40b) correspondingly. At later times we can expand the Green's function (41) in inverse 



powers of t. The first terms in the r.h.s. of (40a) and (40b) are subleasing and we obtain 



V x (x,z,t) 

and by the same token 

V z (x,z,t) i 



1 

P Jo 

■~f 

P Jo 
P Jo 



ds 



ds 



ds 



dr 



1 



dP(r, t 



— OO 

1 



V2-7ra 2 s 



V 4:ira 2 s 



V2na 2 s dr 
[P(R(s),t-s)-P(0,t-s)} 

[P(R( s ),t-s)-P(0,t-s)} 



(43a) 



(43b) 



where R(t) is the boundary beyond which the density of the plasma is below the critical value. We 
observe that V x /V z = Consequently, the azimuthal anisotropy of the hydrodynamic flow is 



V 2 - V 2 



1 - i 

1 2 



V 2 + V 2 l+i 3 (44) 
Since we assumed that the initial conditions and the pressure are isotropic, the azimuthal asym- 



metry (44) is generated exclusively by the magnetic field. 
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IV. SUMMARY 



The structure of the viscous stress tensor in a very strong magnetic field (37) is general, model 
independent. However the precise amount of the azimuthal anisotropy that it generates is of course 
model dependent. We however draw the reader's attention to the fact that analysis of [5] using 
quite different arguments arrived at a very similar estimate. Although a more quantitive numerical 
calculation is certainly required before a final conclusion can be made, it looks very plausible that 
the QGP viscosity is significantly higher than the presently accepted value extracted without taking 
into account the magnetic field effect |14H16j and is perhaps closer to the value calculated using 
the perturbative theory (TTJ HH] . 
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